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ABSTRACT 

Post-translational modifications (such as ubiquitina- 
tion) of clock proteins are critical in maintaining the 
precision and robustness of the evolutionarily con- 
served circadian clock. Ubiquitination of the core 
clock transcription factor BMAL1 (brain and mus- 
cle Arnt-like 1) has recently been reported. How- 
ever, it remains unknown whether BMAL1 ubiquiti- 
nation affects circadian pacemaking and what ubiq- 
uitin ligase(s) is involved. Here, we show that activat- 
ing UBE3A (by expressing viral oncogenes E6/E7) 
disrupts circadian oscillations in mouse embryonic 
fibroblasts, measured using PER2::Luc dynamics, 
and rhythms in endogenous messenger ribonucleic 
acid and protein levels of BMAL1. Over-expression 
of E6/E7 reduced the level of BMAL1, increasing its 
ubiquitination and proteasomal degradation. UBE3A 
could bind to and degrade BMAL1 in a ubiquitin 
ligase-dependent manner. This occurred both in the 
presence and absence of E6/E7. We provide in vitro 
(knockdown/over-expression in mammalian cells) 
and in vivo (genetic manipulation in Drosophila) ev- 
idence for an endogenous role of UBE3A in regulat- 
ing circadian dynamics and rhythmic locomotor be- 
haviour. Together, our data reveal an essential and 
conserved role of UBE3A in the regulation of the 
circadian system in mammals and flies and identify 
a novel mechanistic link between oncogene E6/E7- 
mediated cell transformation and circadian (BMAL1) 
disruption. 



INTRODUCTION 

Circadian clocks drive ~24-h rhythms in gene expression, 
physiology and behaviour and are strikingly well conserved 
in flies and mammals (1,2). At the molecular level, circa- 
dian oscillations rely on transcriptional/translational feed- 
back loops. In mammals, the circadian transcription factors 
BMAL1 (brain and muscle Arnt-like 1) and CLOCK acti- 
vate transcription of two repressors PERIODS (PERs) and 
CRYPTOCHROMES (CRYs) in turn feedback to periodi- 
cally suppress their own transcription (1,3). 

Post-translational modifications modulate protein stabil- 
ity to ensure robust 24-h oscillations. Timely degradation 
of the circadian clock proteins (such as PERs, CRYs and 
REV-ERBa) has been shown to be essential for circadian 
timing (4-8). Degradation by the 26S proteasome requires 
the transfer of ubiquitin moieties to substrate proteins, the 
final step of which requires a substrate specific E3 ubiqui- 
tin protein ligase. It is known that BMAL1 is targeted for 
degradation by ubiquitination (9). However, little is known 
about the ubiquitin ligase(s) targeting BMALland the po- 
tential significance of this to circadian pacemaking. 

Of the ubiquitin ligases, UBE3A (ubiquitin protein lig- 
ase E3A, also known as E6-AP, E6-associated protein) is 
unique in that it can be activated by the E6 oncoprotein 
from high risk human papillomaviruses (HPV16/18), in ad- 
dition to its normal physiological regulators. HPV infec- 
tions underlie several human cancers and are responsible 
for almost all cervical cancers (10,11). During tumorigene- 
sis, the E6/UBE3A complex targets p53 for ubiquitination 
and degradation (12), a major tumour suppressor. However, 
additional E6 targets have been identified that often con- 
tribute to tumour formation (13,14). Core clock proteins 
interact with components of the cell cycle to rhythmically 
control cell-cycle genes, such as Weel, Cyclin Dl, C-myc 
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and p21(Wafl/Cipl ) (15-21), and PERIODS and BMAL1 
have been described as tumour suppressors (16,20,21). To 
date, a direct link between cell transformation and circadian 
clock disruption has not been established. Here, we show 
that E6/E7-mediated cell transformation disrupts molec- 
ular clock rhythms through increasing UBE3A-dependent 
degradation of BMAL1. Further, we provide experimental 
evidence implicating UBE3A as an endogenous key regula- 
tor of the circadian pacemaking in mammals and flies, re- 
vealing the physiological role of UBE3A in the biological 
timing system. 

MATERIALS AND METHODS 

Mice 

PER2::Luc mice were from Prof. J Takahashi (22). All an- 
imal work was carried out in accordance with the 1986 
Home Office Animal Procedures Act (UK) and following 
local ethical review. 

Antibodies and reagents 

The following antibodies were used in this study, BMAL1 
(mouse monoclonal, 23), ubiquitin (Abeam, rabbit poly- 
clonal), REV-ERBa (kind gift from Dr J. Ripperger), 
UBE3A (Cell Signalling, rabbit monoclonal), total GSK30 
and p-GSK3(3 (Cell Signalling, rabbit polyclonal), Myc and 
tubulin (Sigma, monoclonal), Human influenza hemagglu- 
tinin (HA) (Santa Cruz, mouse monoclonal), Histidine tag 
(HIS) (Santa Cruz, mouse monoclonal), V5 (Invitrogen, 
mouse monoclonal), p53 (Merck Millipore, mouse mono- 
clonal) and green fluorescent protein (GFP) (Santa Cruz, 
rabbit polyclonal). MG132, dexamethasone, forskolin and 
cycloheximide were purchased from Sigma. 

Plasmids 

mCLOCK and hBMALl expression plasmids were kind 
gifts from Dr Kazuhiro Yagita (Kyoto Prefectural Univer- 
sity of Medicine, Japan). mAVP::\uc plasmid and luc assay 
were described previously (23,24). Full-length complemen- 
tary deoxyribonucleic acid (cDNA) encoding V5-BMAL1 
was a kind gift from Dr Pat Nolan (MRC Harwell, UK), 
pCW7-HIS-ubiquitin were described previously(25). HA- 
UBE3A (Addgene 8648), HA-UBE3A-C833A (Addgene 
8649), Myc-UBE3A (Addgene 1 1468) and E6/E7 (Addgene 
8641). 

Mouse embryonic fibroblast (MEF) isolation and transfec- 
tion 

MEFs were isolated from E13.5 wild type (WT) PER2::Luc 
mice as described (26). MEFs were transfected with E6/E7 
construct as described before (27) and stable transfec- 
tants confirmed by genotyping and reverse transcriptase- 
polymerase chain reaction (RT-PCR) of E6/E7 fragments 
from genomic DNA and cDNA. NIH3T3 cells were 
purchased from the American Type Culture Collection. 
Cells were grown in complete Dulbecco's modified Eagle's 
medium (DMEM; Invitrogen) supplemented with 10% fe- 
tal bovine serum (FBS) in a 37° C incubator maintained at 
5% C0 2 . 



Soft agar assay 

Agarose was prepared in 6-well plates, with a base layer 
consisting of 1% agar in 1 x DMEM supplemented with 
10% FBS and a top layer consisting of 0.35% agar in 1 x 
DMEM with 10% FBS. Cells were seeded directly into the 
top agar layer at a density of 5 x 10 3 cells per well and al- 
lowed to grow for approximately 3 weeks prior to imaging 
under phase contrast. 

Inimunocytochemistry of p53 

Immunocytochemistry was performed as described (24). 
Briefly, cells on sterile coverslips were fixed with 4% 
paraformaldehyde, then blocked in 1% donkey serum before 
incubation with primary antibody. Cells were then treated 
with an Alexa-Fluor-conjugated secondary antibody mix. 
Fluorescence was visualized using a Zeiss Observer Dl 
AXlO-inverted microscope, photographed using an Axio- 
Cam ICM1 camera and Axio Vision version 4.8.2 software 
(Carl Zeiss). 

Real-time bioluminescence recording and analysis 

Circadian rhythms in bioluminescence from Per2::Luc in 
SW1353 cells or PER2::Luc in MEFs were recorded in real 
time using photomultiplier tube (PMT) devices (28). Con- 
fluent cells were synchronized with 50% horse serum for 2 h, 
Dex (100 nM), Forskolin (FSK) (10 |jlM) or temperature cy- 
cles (36-38. 5°C) as indicated, then transferred to standard 
recording media as described (29). For shRNA knockdown 
in NIH3T3 cells with pGL3-wBwa/7::Luc-Hyg (30), 35-mm 
dishes of cells were transiently transfected by shRNA vec- 
tors. The cells were treated with 100-nM Dex for 2 h and 
then transferred to standard recording media. Biolumines- 
cence was recorded for 1 min at 9-min intervals for 5-10 
days at 37°C in air with Dish Type Luminescencer, Kronos 
(ATTO). Raw data of the rhythms were detrended by 24-h 
moving averages. 

Transfections, immunoblotting and protein stability assay 

Whole cell extracts were prepared in lysis buffer contain- 
ing 1% Triton X-100 as previously described (25). Plas- 
mid transfections for HEK293T or NIH3T3 cells were per- 
formed using Polyfect (Qiagen) according to the manufac- 
turer's instructions. To analyse protein stability, the protein 
synthesis inhibitor cycloheximide (20 |xg/ml) was added to 
cells for the times indicated in the figure legends. Cells were 
then washed and samples analysed over a 6-h period by 
western blot. 

Ubiquitination assay 

Ubiquitination of BMAL1 was examined by histidine 
(His) purification method (25). HEK 293T cells were 
transfected with V5-BMAL1, His-ubiquitin, E6/E7, wild- 
type or mutant derivatives of UBE3A plasmids in varied 
combination as detailed in figure legend. After transfec- 
tion, the cells were left overnight and then treated with 
MG132 for 6 h, followed by solubilizing the cells with ly- 
sis buffer A containing 8-M guanidinium HC1, 100-mM 
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Na2HPO4/NaH2PO4,10-mM Tris-HCl (pH8.0) and 10- 
mM imidazole. Subsequently, Ni-nitrilotriacetic acid purifi- 
cation of the His-ubiquitin conjugates was performed under 
denaturing conditions. 

Co-immunoprecipitation (Co-IP) of Myc-UBE3A, HA- 
UBE3A and V5-BMAL1 

Co-IP of over-expressed protein was performed as de- 
scribed before (25). Briefly, HEK 293T cells in 10-cm dishes 
were transfected with various combinations of tagged plas- 
mids. Overnight after transfection, the cells were solubi- 
lized with lysis buffer (150-mM NaCl, 20-mM Tris-HCl, 
pH7.5, 10% glycerol, 1% Triton X-100, 1-mM phenyl- 
methylsulfonyl fluoride, PhosSTOP and cOmplete ethylene- 
diaminetetraacetic acid-free cocktail protease inhibitors, 
Roche). After centrifugation, 1-ml supernatants were used 
for immunoprecipitation by adding 2 |xg of V5 or HA anti- 
body as appropriate at 4°C overnight. The following day, 
Protein G Dynabeads (Invitrogen) were added and incu- 
bated for 2 h, followed by four washes with lysis buffer in the 
presence of protease inhibitors. Precipitated proteins were 
examined by sodium dodecyl sulphate-polyacrylamide gel 
electrophoresis (SDS-PAGE) and western blot. 

Gene expression analysis 

Ribonucleic acid (RNA) was extracted from cultured 
fibroblasts using the Qiagen RNeasy purification sys- 
tem. cDNA was prepared using the Superscript II re- 
verse transcriptase (Invitrogen) and analysed for gene ex- 
pression using quantitative real-time PCR with TaqMan 
(Applied Biosystems) chemistry (24). Probeset was or- 
dered (pre-validated) from Applied Biosystems (Bmall 
Mm00500226_ml). 

ShRNA and siRNA transfections 

For knockdown of Ube3a, shRNA was designed us- 
ing siDirect (http://sidirect2.rnai.jp/), a web-based soft- 
ware, and the following sequences were used: shUbe3a#l- 
-5'-GGUUU ACUAU GCAAA UGUAG U-3'— and 
shUbe3a#2— 5'-GGAUU AUAUU AUGAC AAUAG A- 
3'. The oligonucleotides to express the shRNA were in- 
serted into the pBS-mU6 vector, as described previously 
(30). Small interfering RNA (siRNA) transfections were 
performed with RNAiMax reverse transfection reagent (In- 
vitrogen). Cells were incubated with 10- or 20-nM SMART- 
pool siRNA duplexes against human Ube3a, or a scrambled 
duplex (Dharmacon) for 24 h before downstream analysis. 
For PMT studies, cells were synchronized with 50% horse 
serum for 2 h before PMT recording. For Western Blotting 
(WB) analysis of knockdown efficiency, transfected cells 
were harvested for protein extraction 48 h post transfection. 

Dvosophila stocks 

All Drosophila stocks were raised on standard cornmeal 
food at 25°C. RNAi line VDRC-45876 was from Vi- 
enna Drosophila RNA interference (RNAi) Center. The 
UAS-ube3a transgenic line was constructed by amplifying 



full-length Dube3a from a cDNA pool derived from whole 
adult flies and cloned into pUAST. Following sequence 
confirmation the UAS-Ube3a was injected into embryos 
for germline transformation. For specific knockdown of 
Dube3a, w;UAS-Gal4;/7^f-Gal4 flies were crossed with 
w;+;\JAS-Ube3a-RNAi to generate w;UAS-Gal4/+;/?#- 
G&l4/UAS-Ube3a-RNAi progeny. For controls each 
line was, instead, crossed to w;+;+ to generate w;UAS- 
Gal4/+;/w//-Gal4/+ and w;+;UAS-C/Z>eJa-RNAi/+. Spe- 
cific over-expression of DUBE3A was achieved by crossing 
u';UAS-Gal4;/>4/-Gal4 and w;+;UAS-Dube3a flies to gen- 
erate H-;UAS-Gal4/+;^/-Gal4/UAS- J D«6eJa. Controls 
(w;\JAS-Gal4/+;pdf-Gal4/+ and w;+;UAS-Dube3a/+) 
were generated by crossing each line to w;+;+. 

Western blotting of dUBE3A in flies 

For the analysis of protein levels, larval brains were dis- 
sected and homogenized in ice-cold lysis buffer [50-mM 
Tris-HCl (pH 7.4), 150-mM NaCl, 1% NP-40, 0.1% SDS 
and 1% proteinase inhibitor (Merck)], and subjected to 
SDS-PAGE according to standard procedures. Western 
blots were probed with anti-UBE3A (8F7, 1:1000, made by 
Yong s lab). The specificity of the antibody was confirmed 
in larval brains from the Dube3a 35 null line. 

Circadian activity assays of flies 

Circadian locomotor activity of flies was measured as de- 
scribed (31). Briefly, male flies, 1-3 days old, were placed 
in glass tubes and entrained to a 12-h light-12-h dark cycle 
(LD) for 5 days and then released into constant dark condi- 
tions (DD) for 7-10 days. Locomotor activity was recorded 
in 30-min bins using DAM2 Drosophila Activity Monitor 
system (Trikinetics, Waltham). Only flies that survived for 
12 days were included in analysis. Data were analysed us- 
ing El Temps software (Antoni Diez Noguera, University 
of Barcelona). Periodogram analysis was done using chi- 
square analysis of LD or DD data, as standard, and cate- 
gorized as arrhythmic if their rhythm strength (variance%) 
score was below the rhythmic threshold. 

Statistical analysis 

Data were evaluated using Student's ?-test or one-way 
ANOVA with Tukey test as indicated. Results are presented 
as mean ± SEM from at least three independent experi- 
ments. Differences were considered significant at the values 
of *P < 0.05, **P < 0.01 and ***P < 0.001. 

RESULTS 

E6/E7-mediated cell transformation disrupts circadian 
rhythm and BMAL1 expression in MEFs 

In vitro, over-expression of HPV16 oncogenes E6 and E7 is 
sufficient to immortalize and transform primary cells (27). 
In a recent attempt to immortalize MEFs isolated from 
PER2::Luc clock reporter mice (22), we used E6/E7 onco- 
genes. Transfection of MEFs with E6/E7 led to success- 
ful integration of viral DNA into the host genome (Sup- 
plementary Figure SI). While mock transfected MEFs no 
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longer proliferated beyond passage 9, E6/E7 expression led 
to indefinite proliferation (beyond passage 30; Figure 1A). 
This was associated with a reduction in p53 protein levels, 
consistent with E6-mediated degradation as previously re- 
ported (Supplementary Figure S2). E6/E7-expressing cells 
were capable of anchorage-independent growth and colony 
formation in soft agar, a hallmark of tumorigenic cell lines 
(Supplementary Figure S3). Thus, stable integration and ex- 
pression of E6/E7 led to immortalization and transforma- 
tion of primary MEFs. Surprisingly, using real-time record- 
ing, the circadian rhythms in PER2::Luc bioluminescence 
were severely dampened in transformed MEFs (Figure IB). 
This was the case in cells synchronized using various well- 
accepted protocols (dexamethasone, horse serum, forskolin 
and 12:12-h temperature cycles; Figure IB and Supple- 
mentary Figure S4). Aberrant rhythms in transformed cells 
were confirmed by studying BMAL1 protein and messenger 
RNA (mRNA) expression following clock synchronization 
(Figure 1C and D). 

BMAL1 availability has been proposed as 'rate-limiting' 
for the positive arm of the clock and is necessary for robust 
rhythmicity (32). In unsynchronized cells average BMAL1 
levels were significantly lower (P < 0.01) in transformed 
cells than in primary cells. This was not accompanied by 
a reduction in steady-state Email mRNA (Figure 2A, P > 
0.05). In addition, the levels of REV-ERBa, a key repres- 
sor of Bmall, did not significantly change (Supplementary 
Figure S5A). These results suggest a post-transcriptional 
mechanism for BMAL1 protein reduction. Indeed, BMAL1 
decay rate was faster in transformed cells than in primary 
cells (Figure 2B) and BMAL1 reduction in transformed 
cells could be rescued by application of the proteasomal 
inhibitor MG132 (Supplementary Figure S5B), suggesting 
the involvement of the proteasomal degradation pathway. 

E6/E7 increases BMAL1 ubiquitination and degradation 
through UBE3A 

BMAL1 is targeted for degradation by ubiquitination 
(9,33,34), thus we sought to address whether E6/E7 expres- 
sion affected this. A Ni-NTA His-ubiquitin pull-down as- 
say was used to precipitate ubiquitinated proteins in HEK 
293T cells. Immunoblotting revealed that E6/E7 expression 
enhanced polyubiquitination of V5-BMAL1 (Figure 2C). 
The final step of ubiquitination requires a substrate-specific 
E3 ubiquitin protein ligase. Due to the well-documented in- 
teraction between E6 and UBE3A (11), we predicted that 
UBE3A may be involved in E6/E7-mediated degradation 
of BMALl. In the presence of E6/E7, over-expression of 
WT UBE3A, but not a dominant negative mutant (C833A) 
that lacks ubiquitin ligase activity, caused a reduction in 
the steady-state level of V5-BMAL1 in HEK 293T cells 
(Figure 2D). These results suggest that the E6/E7-mediated 
BMALl reduction involves the ubiquitin ligase activity of 
UBE3A. 

UBE3A co-purifies with BMALl in a proteasomal pathway- 
dependent manner 

To confirm that BMALl and UBE3A are capable of in- 
teracting with each other in cells, we used Co-IR In the 



presence of E6/E7, V5-BMAL1 was pulled down specif- 
ically by HA-UBE3A and, reciprocally, myc-UBE3A was 
pulled down by V5-BMAL1 (Figure 3A and B). Substan- 
tially more myc-UBE3A was pulled down when the pro- 
teasome was inhibited with MG132, supporting the notion 
that this interaction leads to proteasomal degradation of 
BMALl (Figure 3C). Importantly, V5-BMAL1 could pull 
down myc-UBE3A in the presence and absence of E6/E7 
(Figure 3C). This is distinct from the E6-dependent inter- 
action between UBE3A and p53 (35), and suggests that 
UBE3A acts as an endogenous ubiquitin ligase of BMALl. 
To prove a physiological interaction between endogenous 
UBE3A and BMALl, an endogenous IP was performed 
in unsynchronized SW1353 cells. Endogenous BMALl was 
pulled down specifically by UBE3A using an anti-UBE3A 
antibody (Figure 3D). It is known that GSK3(3 phosphory- 
lates BMALl, promoting its ubiquitination and degrada- 
tion (34). However, we found no difference in the level of 
either total or phosphorylated GSK3|3 (Ser9) between pri- 
mary and transformed cells (Supplementary Figure S6). 

UBE3A regulates clock period and amplitude in mammalian 
cells 

Because BMALl is essential for rhythmic physiology and 
UBE3A interacted with BMALl under physiological con- 
ditions independently of E6/E7, we hypothesized that 
UBE3A-mediated BMALl degradation could be impor- 
tant for endogenous circadian timekeeping. To test this, we 
used shRNAs to knockdown Ube3a in NIH3T3 mouse fi- 
broblasts containing a Bmall ..Luc reporter. Knockdown of 
Ube3a led to a significant reduction in the amplitude of 
Bmall::Luc oscillations (P < 0.01 or 0.001; Figure 4A and 
C and Supplementary Figure S7A). Furthermore, knock- 
down of human Ube3a in a Per2::Luc SW1353 human cell 
line (24) led to unstable clock oscillations, with an initial pe- 
riod lengthening (by up to 4.5 h; P < 0.05) and amplitude 
reduction (-50%; P < 0.05 for 10 nM and P < 0.01 for 20- 
nM siRNA), followed by a loss of (or irregular) rhythmic- 
ity (Figure 4B and C). Consistent with a role of UBE3A in 
regulating BMALl, we observed increased BMALl levels 
following RNAi treatments against Ube3a (Supplementary 
Figure S7B). Ubiquitination of BMALl has been shown 
to be essential for its transactivational activity (33). To es- 
tablish a role of endogenous UBE3A in regulating BMALl 
function, we performed transactivation assay in HEK 293 
cells using a mAVPwhxc promoter reporter (23). The abil- 
ity of BMALl /CLOCK to activate mAVPvXuc was signif- 
icantly reduced by knockdown of Ube3a (Supplementary 
Figure S8), supporting UBE3A as an endogenous regu- 
lator of the core clock complex. Taken together, our re- 
sults demonstrate that loss-of-function of UBE3A in cells 
dampens circadian amplitude and lengthens clock period, 
indicating UBE3A as an endogenous regulator of circadian 
pacemaking in mammalian cells, at least partly through reg- 
ulating the core clock transcription factor BMALl. 

Essential role of dUBE3A in controlling rhythmic fly be- 
haviour 

To determine the importance of UBE3A in driving rhyth- 
mic behaviour, we used Drosophila melanogaster, which 
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Figure 1. Stable integration and expression of E6/E7 immortalized and transformed MEFs and disrupted the circadian clock. (A) Phase contrast mi- 
crograph showing E6/E7 cells were capable of growing beyond passage 10 with morphological change while WT primary cells had become senescent. 
Scale 100 (jim. (B) Representative PER2::Luc bioluminescence rhythms in WT (n = 3) and E6/E7-positive (n = 5) cells following three synchronization 
protocols (dex pulse, forskolin and 50% horse serum). Rhythm dampened rapidly in all cases and period was irregular and significantly longer (28.44 ± 0.77 
h compared to 24.83 ± 0.84 h in primary cells; dex synchronized, P < 0.05). Y-axis value xlOOO = cpm (photon counts per minute). (C) Immunoblotting 
(IB) showing BMAL1 protein rhythm in cells synchronized with dex (100 nM). (D) Bmall gene expression (qPCR) in dex-synchronised cells, n = 3. ***p 
< 0.001; ns, not significant. 
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Figure 2. E6/E7 expression enhanced BMAL1 ubiquitination and proteasomal degradation through UBE3A. (A) Left panel: IB of steady state BMAL1 
in WT and E6/E7 cells. Right panels: integrated density measurement of BMAL1 (n = 3, P < 0.05) and steady state Email gene expression by qPCR (n 
= 3, P > 0.05). **P < 0.01, see underlined. (B) Top panel: BMAL1 protein level following cycloheximide (20 (j,g/ml). Bottom panel: protein degradation 
rate. (C) Ubiquitination assay of V5-BMAL1 in HEK 293T cells by nickel affinity pull-down. Total levels of V5-BMAL1 in the input lysates are shown 
in the panel below. MG132 (5 |xM) was applied in all cases. (D) V5-BMAL1 levels in HEK 293T cells, transfected with E6/E7, co-transfected with either 
HA-UBE3A or a dominant negative mutant (C833A). Representative, n = 3. 



shares high homology with mammals in both the Ube3a 
gene and in molecular clock mechanisms (2,36). We first 
confirmed effective reduction/increase of average dUBE3A 
protein levels in all neurons by driving dUbe3a RNAi 
(VDRC-45876) or dUbe3a-WT (i.e. over-expression; OE) 
pan-neuronally (Figure 5A). Next, we used the pdf>Ga\4 
driver to express dllbe3a RNAi and WT exclusively in the 
LN V neurons that coordinate locomotor behaviour (Fig- 
ure 5B). Both the pdf>dUbe3a-RNAi and pdf>dUBE3A- 
OE lines exhibited normal behaviour under a 12:12-h light- 
dark (LD) cycle (Figure 5C and D). The pdf>dUbe3a- 
RNAi line gradually lost rhythmicity in DD, with 1 1/14 flies 
becoming arrhythmic in DD (Figure 5C and D and Table 1 ). 
Surprisingly, the ;;c/f>dUBE3A-OE exhibited a similar and 



more pronounced arrhythmic phenotype, with 100% of flies 
showing an immediate loss of rhythmicity in DD (Figure 5C 
and D and Table 1). Taken together, our results indicate that 
tight control of dUBE3A levels in central clock neurons of 
flies is critical for the regulation of locomotion rhythms. 

DISCUSSION 

Timely degradation of the clock proteins through the 
ubiquitin-proteasome pathway is critical for the precision 
and robustness of the circadian clock. In contrast to the ex- 
tensive studies on the degradation pathways of clock repres- 
sors (PERs, CRYs and REV-ERBa) (4-7), the degradation 
mechanisms of positive clock proteins (such as BMAL1) 
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Figure 3. UBE3A interacts with BMAL1 both in the presence and absence of E6/E7. (A) Co-IP of HA-UBE3A and V5-BMAL1 in HEK 293T cells. IP 
was performed using an anti-HA antibody and precipitated V5-BMAL1 or HA-UBE3A was detected by IB using an anti-V5 or anti-HA antibody (top and 
middle panels). (B) V5 Co-IP showing that V5-BMAL1 could pull down Myc-UBE3A. MG132 (5 jjlM) was applied in all cases. (C) Co-IP of V5-BMAL1 
and myc-UBE3A in HEK 293T cells in the presence or absence of MG132. IP was performed using an anti-V5 antibody and precipitated myc-UBE3A or 
V5-BMAL1 was detected by IB (top two panels). (D) Endogenous Co-IP between UBE3A and BMAL1 in SW1353 human chondrocytes in the presence 
of MG132. IP was performed using an anti-UBE3A antibody and the precipitated BMAL1 was detected by IB (top two panels). 



are poorly understood. In a recent attempt to immortal- 
ize MEFs isolated from PER2::Luc reporter mice, we used 
E6/E7 oncogenes. Unexpectedly, we observed profound 
disruption of clock gene dynamics in E6/E7-transformed 
cells. Further investigation identified UBE3A as an endoge- 
nous regulator of BMAL1 ubiquitination and degradation, 
controlling molecular and behavioural rhythmicity in mam- 
malian cells and in flies. 

Consistent with the near ubiquitous expression of circa- 
dian clocks, UBE3A is widely expressed in multiple mouse 
tissues, including the central pacemaker superachiasmatic 
nucleus (SCN) (37 and Supplementary Figure S9). Both 
over-expression and knockdown of UBE3A in fly central 
clock neurons resulted in a pronounced circadian pheno- 
type, indicating that precise control of its activity is vi- 



tal for circadian pacemaking across divergent species. It 
has been reported that flies carrying a global null muta- 
tion of dUbe3a demonstrate circadian deficits in locomo- 
tor rhythms (36). The clock-neuron-specific manipulations 
presented here provide unequivocal evidence that dUBE3A 
is an important regulator of the Drosophila central clock 
and suggests a conserved role of UBE3A in circadian time- 
keeping in flies and mammals. It is possible that dUBE3A 
may target CYCLE, the fly orthologue of BMAL1. How- 
ever, the lack of a commercially available antibody for CY- 
CLE prevented us from testing this possibility. Of note is 
that UBE3A plays a causal role in a human neuro-genetic 
disease, Angelman syndrome (38,39). This neurological dis- 
order is characterized by severe intellectual and develop- 
mental disability, seizures and sleep disturbance. It is un- 
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Figure 4. Knockdown of endogenous UBE3A disrupts circadian rhythms in mammalian cells. (A) Representative raw (left) and detrended (right) biolumi- 
nescence traces of Bma/l::Luc NIH3T3 cells treated with shRNA #1 (top) and #2 (bottom) against UbeSa. (B) Representative traces of Per2:'A\iC SW1353 
cells treated with human Ubeia SiRNA constructs at 10 (mid-grey) and 20 nM (black) doses. Scrambled SiRNA was used as control. (C) Left panels: 
circadian amplitude analysis (n = 3, shRNA #1 P < 0.01; shRNA #2 P < 0.001; RNAi 10 nM, P < 0.05; RNAi 20 nM, P < 0.01). Right panel: circadian 
period analysis (n = 3, t-tests: RNAi 10 nM, P < 0.05; RNAi 20 nM, P < 0.05). Note: the rapid dampening of circadian oscillation in shRNA-treated 
NIH3T3 cells prevented proper evaluation of periodicity. *P < 0.05; **P < 0.01; ***p < 0.001. 
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Figure 5. Drosophila dUBE3A is necessary for circadian rhythmic behaviour in central clock neurons. (A) IB of total protein from L3 larvae brain, 
showing complete lack of protein in the dUBE3A-mi\\ (DUBE3A 35 ), efficiency of dUBE3A-RNAi knockdown (RNAi) and over-expression (OE) by the 
pan-neuronal driver e/av-Gal4. Two independent blots showed similar results. (B) Expression of GFP driven by the pdf-Ga\4 driver, showing restriction to 
the ventrolateral neuron (LN V ) subgroups of central clock cells in the fly brain that coordinate daily sleep-wake cycles. (C) Upper panels: representative 
double-plotted actograms for locomotor activity data in flies with Pigment dispersing factor (PDF) neuron-specific dUBE3A RNAi or dUBE3A over- 
expression, along with relevant controls. LD: white background; DD: grey background. Light and dark bars at the top indicate day (ZTO-12) and night 
(ZT12-24) during LD, or, subjective day (CT0-12)/subjective night (CT12-24) for DD, respectively. Lower panels: representative periodograms during 
LD and DD for each genotype. (D) Relative percentage of rhythmic flies for each genotype; n = 10-33. 



Table 1. The circadian phenotype of transgenic fly lines expressing dUbe3ci RNAi (light grey) and over-expressing dUBE3A (dark grey) 



Mean ± SEM 





Arrhythmic 
in LD 


Arrhythmic 
in DD 


DD period 


DD % 
variance 


n 


pdf-gal4 
driver control 


0 


0 


23.65 ± 0.32 


36.64 ±1.56 


28 


RNAi 

responder 

control 


0 


0 


23.96 ± 0.92 


39.88 ± 3.63 


32 


pdf-gal4 x 
UAS-RNAi 


0 


28 


22.53 ±1.35 
(based 
on 5 flies) 


23.46 ±2.08*** 


33 


UAS- 
dUBE3A 
Responder 
Control 




0 


24.49 ± 0.29 


28.41 ±1.61 


10 


pdf-gal4 x 

UAS- 

dUBE3A 


0 


11 


n/a 


14.59 ±1.03*** 


11 



***P < 0.001. 

Period, h; % variance denotes rhythm strength. 
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known whether the dis-regulation of the UBE3A-mediated 
circadian rhythms might contribute to some of the symp- 
toms of this human disorder. 

BMAL1 ubiquitination has been suggested as a possible 
global mechanism of clock entrainment (9,33,34). However, 
our study is the first to report a ubiquitin ligase targeting 
BMAL1. Interestingly, a recent study revealed UBE3A as 
a light responsive protein in the mouse SCN (36). It re- 
mains to be seen whether UBE3A is involved in light or 
food-resetting of the circadian clocks. The core clock pro- 
teins PERIODs (and BMAL1 more recently) have been 
described as tumour suppressors, rhythmically controlling 
cell-cycle genes (16,20,21). For instance, downregulation 
of BMALl promotes tumour growth in vivo and acceler- 
ates cancer cell invasion in vitro, partially through its in- 
terplay with the tumour suppressor p53 pathway. In a re- 
cent genome-wide RNAi screen, downregulation of Bmall 
allowed cells to bypass p53-dependent cell-cycle arrest (20). 
Our data show that expression of HPV16 E6/E7 oncogenes 
causes increased degradation of BMAL1 and disrupts circa- 
dian oscillations. Although extensive further work is needed 
to draw a definitive link, it is tempting to speculate that 
this disruption might affect the expression of downstream 
clock-controlled genes regulating the cell cycle (such as Cy- 
clin Dl and Weel), contributing to tumorigenesis in HPV- 
related cancers. 

Together, our results indicate UBE3A, an E3 ubiquitin 
ligase previously implicated in HPV-tumorigenesis and An- 
gelman Syndrome, as a conserved and integral part of the 
circadian clock feedback loop. These findings not only re- 
veal a de novo key regulator for the circadian transcription 
factor BMAL1 but also offer a new mechanistic link be- 
tween cell transformation, intrinsic clock disruption and 
cancer. 
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